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Dependence of the energy of ions accelerated during interaction of the laser pulse obliquelly 
incident on the thin foil target on the laser polarization is studied experimentally and theoretically. 

We found that the ion energy being maximal for the p-polarization gradually decreases when the 
pulse becomes s-polarized. The experimentally found dependences of the ion energy are explained 
by invoking the anomalous electron heating which results in high electrostatic potential formation 
at the target surface. Anomalous heating of electrons beyond the energy of quiver motion in the 
laser field is described within the framework of theoretical model of driven oscillator with a step-like 
nonlinearity. We have demonstrated that the electron anomalous heating can be realized in two 
regimes: nonlinear resonance and stochastic heating, depending on the extent of stochasticity. We 
have found the accelerated ion energy scaling determined by the laser intensity, pulse duration, 
polarization angle and incident angle. 


I. INTRODUCTION 

Discovery of laser-driven ion acceleration [lj has 
opened a broad field of potential applications involving 
fast ignition ii and hadron therapy ■ Most inten¬ 
sively investigated has been the mechanism of ion acceler¬ 
ation with charge separation (CS) fields generated on the 
thin-foil target surfaces, including target normal sheath 
acceleration (TNSA) jl|. Nowadays, increasing atten¬ 
tion is also paid for different mechanisms, however, the 
TNSA mechanism still remains crucial to understanding 
physics of ion acceleration at the moderate laser intensity 
regimes. 

Generation of ion accelerating CS field is governed by 
the absorption mechanism of laser energy into electrons. 
In the case of interactions between obliquely incident 
laser and overdense plasma, one of the most predomi¬ 
nant absorbed mechanism is the Brunei effect [7] . In this 
model, electrons on the boundary are directly accelerated 
by the electric field of the laser, when the electrons are 
accelerated to the normal direction of the surface with a 
momentum p max = 2m e cag at maximum, where c is the 
speed of light in vacuum, m e and e are the electron mass 
and charge; ao = eEi /ra e wc is a dimensionless amplitude 
of the laser radiation, with Ei and lo being the laser elec¬ 
tric field amplitude and frequency. It can be expressed in 
terms of the laser intensity / and wavelength A = 2i tc/w 

as o 0 = O^y^lO 18 W/cm 2 )A 2 (/mr 2 ). 

To explain the acceleration beyond p m ax , nonlinear ef¬ 
fect should be implemented into the analysis. D’yachenko 
and Imshennik [8( numerically revealed that electrons 
quivered by the laser field are kicked-out by the potential 
boundary of the plasma and anomalously gain momen¬ 
tum depending on time (see also 0E2).' Paradkar et al. 
P] and Krasheninnikov [10] theoretically showed stochas¬ 
tic heating of electrons by assuming a model dependence 


on coordinate of the electric field corresponding to the 
field in the vicinity of a thin-foil tar gets . Using particle- 
in-ccll simulations, Taguchi et al. 13] found that the 
heating of electrons transiting through a small size tar¬ 
get can be higher than that by the electron quivering 
and called this phenomenon as “nonlinear resonance ab¬ 
sorption.” Robinson et al. [l4| proposed non-wake-field 
acceleration beyond the momentum limit is gas plasma. 

In this manuscript, we develop the mechanism of ion 
acceleration involving nonlinear effect on the electron en¬ 
ergization in order to describe the accelerated proton en¬ 
ergy during the laser interaction with the thin foil tar¬ 
get. By considering a model of driven oscillator with a 
step-like nonlinearity, the growth of chaotic behaviour of 
electrons is revealed. Here, we call the mechanism the 
Nonlinear Vacuum Heating. As a result, the experimen¬ 
tal results on the proton energy and its dependency on 
the laser polarization are successfully predicted. 


II. EXPERIMENTS 

The experiment was performed using JLITE-X laser 
at JAEA, which delivers ultrashort (r = 40 fs, FWHM) 
linearly-polarized pulses having a central wavelength of 
A = 800 nm. The laser contrast has been raised to 10 9 
by the benefit of Insertable Pulse Cleaning Module [lf| . 
The beam is focused to a spot of 8 ^m (FWHM) in di¬ 
ameter using an / = 160 mm off-axis parabolic mirror 
(OAP) with an incidence angle of 45° on silicon-nitride 
foils. The polarization angle is varied from p (0 = 0°) to 
s (0 = 90°) by using a half-waveplate located at the en¬ 
trance of the OAP. The areal peak intensity obtained is 
I = 3.0 x 10 18 W/cm 2 , which corresponds to the dimen¬ 
sionless amplitude ao = 1.2. The pedestal intensity level 
in sub-ns range is around 10 9 W/cm 2 that is adequate 
to suppress the preplasma formation on the ultrathin foil 
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P P° ■ Polarization angle [degree] s pol. 

FIG. 1. Maximum proton energies (circles) observed for 500 
nm thick targets as a function of the polarization angle 8. 
When the laser polarization is close to s-polarization (8 > 
80°), the proton energy becomes below than the detection 
limit. Analytically predicted proton energy £ p is shown with 
a solid line. The inset shows a typical energy spectrum for 
8 = 0 . 


surface. Accelerated protons are measured by Time-of- 
Flight (TOF) online detector m located at the normal 
direction of the target rear surface. It is confirmed us¬ 
ing ion-track (CR-39) detection that the protons (> 1.2 
MeV) has a 15° angular divergence distributed along the 
normal of the target surface [16|. 


In the experiment we analyze the dependence of the 
energy of accelerated ions on the thin foil target on the 
laser polarization. Figure [0 shows the maximum energy 
of protons (circles) observed for 500-nm-thick targets as 
a function of the polarization angle 9. Here, the pro¬ 
ton energy is drastically decreasing as the polarization is 
changed from p to s mode. Similar dependency was also 
observed [13 for a thicker (13 /im) target under high- 
contrast laser conditions. On the other hand, there is a 
clear difference from the case of low-contrast laser when 
the proton energy for s-polarization was only by 30% 
lower [13 than that for p-polarization. As we see, the 
obvious dependency on the polarization angle is charac¬ 
teristic of the interaction between steep-gradient plasma 
surfaces and obliquely-incident laser pulses. To under¬ 
stand the physics underlying the phenomenon, we discuss 
the motion of electrons driven by a strong EM wave in 
the vicinity of the plasma-vacuum interface. 


(a) y (b) y' (c) E 


E' 



FIG. 2. Electromagnetic wave obliquelly incident on the thin 
foil target, (a) The wave in the laboratory frame of reference 
K. (b) The wave in the boosted frame of reference K'. For 
the p-polarization wave, the electric field vector is in the x-y 
plane. In the boosted frame of reference, the plasma electrons 
move along the target with a velocity — V = — csin tj> in the y' 
direction, (c) A charge separation electric field on the x axis. 


III. OBLIQUE INCIDENCE OF THE LASER 
RADIATION ON A THIN FOIL TARGET 

A. Electron quiver motion 

Let us consider a plane electromagnetic (EM) wave 
interacting with the thin foil located in the plane x = 0. 
The laser the angle of incidence equals <^> (see Fig. [2] 
(a)). The EM wave is linearly polarized and its electric 
field E has a polarization angle 8 with respect to the 
x-y plane: 9 = 0 for p-polarization and 8 = 90° for s- 
polarization. In order to simplify the consideration, we 
perform a Lorentz transformation from the laboratory 
frame K to the frame of reference K' moving with the 
velocity V = c sin (j) along the target surface (y direction). 
In the K' frame, the EM wave is incident normally onto 
the surface [U (22|, and the gamma factor associated 
with the velocity V is given by 

7M y/1 - U 2 /c 2 COS0 ^ 

Fourier component for a monochromatic plane wave in 
vacuum [l9| are given in K' frame by equations 

VxE' = i— B' (2) 

c 

V x B' = -i— E', (3) 

c 

where the wave frequency equals u / = w/tm = uj cos (/). 
Using V • E' = 0, the Eqs. © and (0 can be written as 
follows 

w ' 2 

V 2 E' + -5-E' = 0. (4) 

c 2 

Using Eqs. ([2] SI) we obtain the general solution of the 
Maxwell equations, 

E' = E 0 e"' (4 ' _3: ' /c) + E r e i “' (t ' +a: ' /c) 


( 5 ) 
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with the electric field amplitudes of incident Eo and re¬ 
flected E r waves in the K' frame. 

When the EM wave interacts with the plasma having 
a steep gradient satisfying the condition r E /L 1, we 
can use the Leontovich boundary conditions HU for the 
tangential components of the electric and magnetic fields 
E t and B t : 

E t = CB t x n, (6) 

where n is the unit vector along the inward normal to the 
plasma surface and £ is the surface impedance operator. 
Here r E is the electron quiver radius and L is the scale- 
length of the plasma inhomogeneity. Assuming that the 
plasma surface can be treated as an ideal conductor (£ = 
0), we obtain the boundary condition 


B. Nonlinear mechanism of anomalous electron 
heating 

Within the framework of the Nonlinear Vacuum Heat¬ 
ing mechanism, the electrons quivered by the laser field 
are kicked-out by the steep-like electric field at the plasma 
boundary. As the result, the EM field energy is irre¬ 
versibly converted into the electron kinetic energy, thus 
enhancing the electron heating. To describe the nonlin¬ 
ear motion of electrons, we assume that a step-like CS 
electric field is formed near the ion layer having the pos¬ 
itive electric charge. The equations of electron motion 
can be written in the form 

p + £sign(x) = aocost, (14) 


K'=o = E t = 0. (7) 

Then, the electric field is written as 

E' = E 0 [e iu;,< ' t '- x 'S c ' > — e iu '( t '+ x ’/ c)j_ (g) 

The magnetic field is obtained using Eqs. m and m as 
follows: 

B' = E 0 [e ia/( *'- a:,/c) + g^'d'+x'/c)^ ( 9 ) 


X = 


V 

(i+p 2 ) 1/2 ' 


(15) 


Here the sign function is equal to sign(x) = 1 for x > 0 
and sign(x) = —1 for x < 0, and a dot denotes a differen¬ 
tiation with respect to time, t. The electron momentum 
is normalized on m e c , the time is measured in units w -1 , 
and the coordinate is normalized on c/ui. 

The parameter 


e p = 2'Knoe 2 l/m e uJC, 


(16) 


Therefore, the real part of the magnetic field at the sur¬ 
face (x 1 = 0) is given by 

B^./_ 0 = (0, — 2E 0 sin# cos a/i', 2Eo cos 6 cos u'f'). (10) 

In the boosted frame of reference, K\ the acceleration 
of electrons in the x' direction is driven by the Lorentz 
force equal to — (e/c)v x B, where the magnetic field is 
taken at the x = 0 plane, B = B( c , =0 , and the velocity 
of the electrons is equal to v= (0, —V, 0). As the result 
the equation of motion for electron can be written as 

% = < u > 

It yields for the ^'-component of the electron momentum 

-E?L = 2m e Vaou/cos0cosu/t'. (12) 

at 

Hence, the x'-component of the electron momentum de¬ 
pends on time as 

p x > = 2m e Vao cos# sin u/i 7 , (13) 

where is the dimensionless amplitude of the laser and 
V = c sin <f>. We see that the electron momentum p x i 
is below 2m e caa , which corresponds to the momentum 
limit for electrons accelerated by obliquely incident laser 
via the Brunei effect [jj. 


introduced in Ref. [22J, is proportional to the charge 
separation electric field, E 0 = 2Trenol (see Fig. [2] (c)), 
normalized on m e a;c/e, for a thin foil with the electric 
charge surface density equal to en^l. 

The formulated model of a driven oscillator with a 
step-like nonlinearity (see also HU) helps to determine 
and elucidate the conditions when the stochastic vac¬ 
uum heating takes place. The 3-dimensional model used 
bears the key features of the problem addressed in Refs. 
[9;, 10, 24, 25J, but in a significantly simpler mathematical 
approximation, because a full description of the charged 
particle interaction with electromagnetic waves implies a 
consideration of a 7-dimensional dynamical system. 

While the electron trajectory crosses the x = 0 plane, 
the sign of the static electric field changes abruptly. This 
can be considered as a “collision”, during which the os¬ 
cillating driver electric field produces the work changing 
the electron energy. Depending on the phase, when the 
“collision” happens, the work can be either positive or 
negative, resulting in the electron energy increasing or 
decreasing. 

The enhancement of the energy transfer from the 
driver electric field to the electron is expected to occur 
when a resonant condition (for details see Ref HU) takes 
place: 

a 0 = tte/2, (17) 

which can be equivalent to the relationship of E 0 « £). 
The expression for the momentum maximum is found in 
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Ref. (23| for two different nonlinear regimes: nonlinear 
resonance and stochastic heating. 

In nonlinear resonance regime, electrons gain a mo¬ 
mentum A p « a 0 during the passage of the plane x = 0. 
Then, the momentum grows proportionally to the square 
root of time, written in 


Pm ~ (a 0 et/2) 1/2 = a 0 (t/n) 1/2 (18) 


p ( a ) 



for the regime of nonlinear resonance. In the case of 
stochastic heating, the momentum dependence on time 
has a character of diffusion with the diffusion coefficient 


2ApAp Oq£ 
T 2p m 


(19) 


This yields for the time dependence of the maximum elec¬ 
tron momentum 


p m « (ag£f/2) 1/3 = a 0 {t/n) 1/3 , (20) 

noticed in 

To analyze in more details the properties of nonlinear 
oscillations we have integrated Eqs. m usd numerically. 
For regularization of the singularity in the left hand side 
of Eq. (fl4l) we use instead the sign function sign(a;) the 
function Tanh(x/Z) with the ion layer width l (see Fig. 
[2] (c)) substantially small compared to the particle dis¬ 
placement: l < Ax. Figures [3] (a,b) show the result of nu¬ 
merical calculation of the momentum and the coordinate 
of electrons as functions of time obtained for the condi¬ 
tion close to our experiment. We can see clearly that the 
oscillations are nonlinear with the momentum abruptly 
droping down and rising up again. In a time duration of 
the laser pulse, 40 fs (t = 95 in Fig. [3j , the momentum 
reaches 8 m e c, which is a few times higher than the quiver 
momentum corresponding to the Brunei effect. The os¬ 
cillation amplitude is limited both in the momentum and 
coordinate directions, shown in Fig. 13c). The trajectory 
tightly fills the domain in the phase plane demonstrat¬ 
ing ergodic property of the system under consideration. 
In the Fig. [3] (d) we plot a Poincare section having a 
form of finite thickness web, which is broadened due to 
stochastic nature of the particle motion. The property of 
the particle to migrate on the phase plane along the web 
is typical for the minimal chaos regimes {HI. In other 
words, chaotic behavior has begun to emerge, but not 
evolves completely under these conditions. 



FIG. 3. Driven ocillations obtained for ao = 2 and e v = 1.35. 
Normarized momentum (a) and coordinate (b) of electrons 
as vs time, (c) The phase plane ( x,p ). (d) The Poincare 
section showing the particle positions on the phase plane at 
the discrete time with the time step equal to the period of the 
driving force, 27 t. 


C. Ion acceleration 

Nonlinear Vacuum Heating conjecture is equivalent to 
the assumption that the energy conversion from the EM 
field to electron energy depends on time according to 
the nonlinear mechanism introduced above. The average 
square of aZ-component of the electron momentum p nv h,x' 
scales with the laser pulse duration which is proportional 


to the number of wave periods N for the two different 
regimes of nonlinear vacuum heating as: 

Np 2 x , (nonlinear resonance) (21) 

N 2 / 3 p 2 , (stochastic heating) (22) 

To estimate the typical energy gain of fast protons, we 
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use the step-like CS field equal to Eo tanh(a;/i) in the 
interval 


— 1/2 — Ax < x < 1/2 + Ax, (23) 

where Ax is the electron displacement and l is the thick¬ 
ness of positively charged layer of the foil target (Fig. 
(21(c)). When the target thickness is thin enough to sat¬ 
isfy the condition l ~ Aa;, the electric field E$ can be 
expressed as 


Eo = 2imoel = 27moeAa: (24) 

with the ion charge density equal to no . Assuming that 
protons can be treated as test particles initially located 
on the surface of the layer, the proton energy gain £ p can 
be estimated to be of the order of the electron kinetic 
energy, £ e \ 


An 2 

£ p = eEoAx = 27rnoe 2 Aa; 2 ~ -— = £ e . (25) 

In the boosted frame of reference K', the kinetic energy 
can be found to be given by 

£' v = £' e = £' total - m e c 2 (26) 

S'totai = 7^ j> y/m 2 e c 4 + P 2 v c 2 + p 2 nvhtX ,c 2 dt, (27) 

where 

Pv = 7 MWeV = m e c tan cj>. (28) 

From Eqs. (fT3l) . (I2T1) and (l22l) . follows that 

Pnvh,x' = N J p x i = 2N 3 m e aoC sin <j> cos 0 sin uj't' (29) 

with the index j equal to j = 1/2 for the nonlinear res¬ 
onance regime and equal to j = 1/3 in the stochastic 
heating regime. The kinetic energy gain in the labora¬ 
tory frame of reference, £ p can be found by performing 
the inverse Lorentz transformation from the boosted to 
laboratory frame of reference, 

£ P = lM{£'totai ~ PvV) - m e c 2 . (30) 


This the equation for the proton energy £ p . Its analytical 
solution yields 


f — 

„„ 2 


TTtpC 2 f 1 


COS^ 0 7T 


E(A) + E 


A 


1 T A 


\/l + A 


- 1 


where 


and 


A = ag N 23 cos 2 6 sin 2 2cj) 


W 2 ---— 

E{z)= / vl-z 2 sin 2 0d0 

Jo 


(31) 

(32) 

(33) 


is the complete elliptic integral of the second kind 



FIG. 4. Maximum proton energy observed for different val¬ 
ues of the the target thickness: 0.1-0.5 /im thick silicon nitride 
and 5-12.5 /jrn thick polyimide. Bars correspond to a double 
standard deviation (2 a) for each point. The value of the pro¬ 
ton energy £ p shown by a solid line is found from Eq. mb¬ 


it should be emphasized that the present model deter¬ 
mines the maximum energy of protons as a function of the 
laser parameters: dimensionless amplitude a 0 , number of 
wave period N, angles of polarization 6 and incidence (j). 
In Fig. 1, the analytically value of proton energy £ p pre¬ 
dicted with nonlinear resonance mechanism (see Ea. (131l) 
with j = 1/2) is shown as a function of 6. Here, we are 
using parameters employed in the experiment: ao = 1.2, 
<fi = 45°, and N = 15, which corresponds to the number 
of wave periods included in the FWHM duration (r = 40 
fs) of the incident laser. On the other hand, the pro¬ 
ton energy is also reproduced by stochastic heating (See 
Ea. (l3T1) with j = 1/3) using (ao ,4>,N) = (1.7,45°, 15), 
where the ao value is slightly higher than the experimen¬ 
tal one. Although it may be possible to assume that the 
enhancement of a 0 from 1.2 to 1.7 is caused by a self- 
focusing of laser in preformed plasma, it is natural to 
consider that the electron energization in the present ex¬ 
perimental condition is predominantly governed by non¬ 
linear resonance mechanism, not by stochastic heating. 
This aspect is consistent with the fact revealed by nu¬ 
merical analysis (Fig. 3) that the process is in the inter¬ 
mediate stage between nonlinear resonance and stochas¬ 
tic heating and that chaotic behavior of energization has 
just begun to grow. Substantially, the chaotic regime re¬ 
quires a long time to develop as an asymptotic stage after 
the nonlinear resonance. 

As is shown in Eg. EH) , the present model is valid when 
the target thickness satisfies l ~ Ax. The displacement 
of electrons Aa: is typically assumed to be equal to the 
amplitude of electron oscillation = eEi/m e u> 2 . Hence, 
the present model is applicable when the target thickness 
satisfies l ~ te ■ To demonstrate applicability of the crite¬ 
rion above, we observed the target thickness dependency 
of proton energy, shown in Fig. (J) 

Fig. [I] shows the dependence of the proton energy 
on the target thickness observed with (ao,0,<f>,N) = 
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(a) nonlinear resonance 
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(b) stochastic heating 



FIG. 5. Isocontours of the proton energy £ p measured in a 
unit of MeV as in the plane of the laser amplitude a o and 
normalized duration N in two different regimes of nonlinear 
vacuum heating: (a)nonlinear resonance and (b) stochastic 
heating. The laser incidence and polarization is fixed on tj> = 
60° and 9 = 0. 

(1.2,0,45°, 15). The maximum proton energy is close 
to the value of analytical prediction (2.18 MeV) given 
by Eq. (HU) when the target thickness is ranging around 
t\e = ao\/2ir = 0.15 pm. For thicker targets (l > 5pm), 
proton energy decreases with the increase of the target 
thickness l. This can be attributed to the neutralization 
of ion density no by the returning current, whose effect 
is subtantially stronger in the case of thicker targets. 


D. Ion energy scaling 

Figure [5] shows the acquired proton energy according 
to Eq. m, under the nonlinear resonance (j = 1/2) (a) 
and stochastic heating (j = 1/3) (b). Isocontours of the 
proton energy £ p measured in a unit of MeV as in the 
plane of the laser amplitude a o and normalized duration 
N. In this figure, the proton energy is constant along a 


contour line OgV 2 ^ =constant. It is interesting to note 
that the nonlinear resonance provides the energy scaling 
on clqN , which is proportional to the laser energy fluence. 

As is discussed above, the stochastic heating requires 
a substantially long time to develop being an asymptotic 
stage after the nonlinear resonance. Therefore, when the 
laser pulse duration is relatively short, typically r ~ 50 
fs ( N ~ 20 for A = 800 nm) as used in our experi¬ 
ment, the scaling perhaps correspond to the nonlinear 
resonance regime. According to this hypothesis, the en¬ 
ergy gain up to 200 MeV (this is one of the milestones 
toward clinical applications of the laser accelerated ions 
0), can be achieved for the laser radiation with the am¬ 
plitude ao = 42 (see Fig. Oja)). This laser amplitude 
corresponds to the intensity of I = 4 x 10 21 W/crn 2 and 
requires the target with optimal thickness of the order 
of l ~ Te = 5.4 pm. On the other hand side, a longer 
pulse duration, e. g. when r = 500 fs for laser wave¬ 
length o f A = 1.06 pm, N = 142, the proton energy 
equal to 200 MeV can be obtained with the laser ampli¬ 
tude equal to a o = 16 corresponding to the laser intensity 
I = 3 x 10 20 W/cm 2 and the target thickness of the order 
of l — 2.6pm, provided the stochastic heating mechanism 
does not occur. The laser parameters written above have 
been achieved in several laser facilities nowadays. How¬ 
ever, for a longer pulse duration (i.e. for r = 500 fs, 
A = 1.06 pm N = 142) in the stochastic heating regime 
the 200-MeV proton acceleration requires the laser ampli¬ 
tude to be equal to a o = 34 (see Fig. [51(b)) corresponding 
to the intensity I = 1.6 x 10 21 W/cm 2 . 

IV. CONCLUSION 

In conclusion, We have studied the ion acceleration 
during interaction of the short pulse laser obliquelly in¬ 
cident on the thin foil target. The dependence of the ion 
energy of on the laser polarization has been studied stud¬ 
ied experimentally and theoretically. We found that the 
ion energy being maximal for the p-polarization gradu¬ 
ally decreases when the pulse becomes s-polarized. The 
experimentally found dependences of the ion energy are 
explained by invoking the anomalous electron heating, 
which results in efficient ion acceleration at the target 
surface. 

We showed that the proton energy can be substantially 
enhanced under the conditions of nonlinear resonance ex¬ 
perienced by the electrons circulating around the thin foil 
target irradiated by the electromagnetic wave. 
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